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A disorder parameter detecting dual superconductivity of the vacuum is measured across the chiral phase 
transition in full QCD with two flavours of dynamical staggered fermions. The observed behaviour is similar to 
the quenched case. 



1. Introduction 

Confinement is experimentally a well estab- 
lished property of strong interactions. In 
the quenched approximation a phase transition 
driven by temperature from a confined to a de- 
confined phase has been observed by lattice simu- 
lations, an order parameter being the v.e.v. of the 
Polyakov line, which is zero in the confined (low 
temperature) phase and nonzero in the decon- 
fined (high temperature) phase, thus signalling 
the spontaneuous breaking of the Z3 centre sym- 
metry. The situation is different in full QCD, 
where the presence of dynamical quarks makes 
the centre symmetry explicitely broken, so 
that the v.e.v. of the Polyakov line ceases to 
be an order parameter. One speaks instead of 
a phase transition from a spontaneously broken 
to a restored chiral symmetry, the chiral conden- 
sate being the order parameter. Anyway, at the 
chiral phase transition a peak in the Polyakov line 
susceptibility is still observed, indicating a sort of 
deconfining crossover. However up to now the re- 
lation between chiral symmetry and confinement 
is still not well understood. 

On the other hand, confinement in the 
quenched theory is now well understood in terms 
of an order/disorder phase transition 0,0], with 
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the QCD vacuum behaving like a dual super- 
conductor at low temperatures and undergoing a 
phase transition to a normal conducting state at 
the critical temperature Tc. A disorder parameter 
which is the v.e.v. of a magnetically charged oper- 
ator has been developed and used to demonstrate 
dual superconductivity of the vacuum for U(l), 
SU(2) and SU(3) Yang-Mills theory Con- 
trary to the v.e.v. of the Polyakov line, the defini- 
tion of the disorder parameter detecting dual su- 
perconductivity can be extended to full QCD. It is 
then possible to test whether there is any transi- 
tion from a dual superconducting to an ordinary 
vacuum also in presence of dynamical fermions 
and how it is related to the chiral phase transi- 
tion. This is the aim of our study and we will 
present some preliminary results in this paper. 

2. The disorder parameter in full QCD 

The definition of the disorder parameter in pure 
gauge theories has been discussed in , so we 
will only recall here the main points. It is the 
v.e.v. of a gauge invariant operator, /i, which 
carries magnetic charge. The U{\) magnetic sub- 
group is selected by the so called abelian pro- 
jection. A gauge fixing is performed, usually 
by diagonalizing a local operator in the adjoint 
representation, which leaves a U(l)xU(l) gauge 
freedom (for SU(3)), corresponding to two pos- 
sible definitions of the magnetic charge. It has 
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been proved in that the dual superconductivity 
mechanism for confinement works independently 
of the abelian projection chosen, as already sug- 
gested in 1^. 

On the lattice ^ can be written as [^j-Q 



exp^J2ReTT(Uo^i£,to)-Uo^{x,to)) , (1) 



where the sum in the exponential is extended only 
to the temporal plaquettes C/oi in the time slice 
t = to where the magnetic monopole is created, 
and Uoi{x,tQ) indicates a plaquette modified by 
the insertion of a monopole field, defined as an 
abelian field in the abelian projected gauge. The 
disorder parameter is then defined as 



/ (VU) ^e-^c _ / (VU) e-^"-' 



im)' 



(2) 



where Sq is the standard Wilson action and Sm 
is a modified action obtained by inserting the 
monopole field in the temporal plaquettes on time 
slice to- It can be shown by a change of 
integration variables, that the insertion of fi in 
the functional integral corresponds to creating a 
monopole at to which then propagates at all times 
t > to- For this reason, periodic boundary con- 
ditions in time direction are not consistent with 
the definition of (/i), and C*-periodic boundary 
conditions have to be used instead: 



U^,ix,t^Nt) = U;{x.,t^O) 



(3) 



Nt being the temporal extension of the lattice and 
U* being the complex conjugated of U^. 

A lattice determination of {fi) is difficult, since 
it is the average of a quantity which fiuctuates 
exponentially with the square root of the physical 
volume. It is much more convenient to study 

d 

P=-Qp log(M) = i^G)sG - {Sm)sm , (4) 

where {■)s indicates the v.e.v. obtained using 
the action S to weight configurations. From the 
knowledge of p the relevant physical information 
on (/i) can be extracted. The result in pure gauge 
theory is that (fi) ^ in the confined, low tem- 
perature phase, while it goes to zero at Tc, with 



critical indices which can be extracted from the 
finite size scaling behaviour of p. 

In full QCD the definitions of (p) and p are 
modified in a very simple way: 



where Sp is the fermionic action, and 

P = {Sg)sg+Sf - {Sm)sm+Sf ■ 



(5) 



(6) 



Unlike the centre symmetry, the U(l) mag- 
netic symmetry defined after abelian projection 
is still a good symmetry also in presence of dy- 
namical fermions. Therefore (p) can be a correct 
disorder parameter for the dual superconductiv- 
ity transition also in full QCD. 

3. Algorithm implementation and simula- 
tion details 

We have used two flavours of staggered 
fermions and the R version of the HMC algorithm 
for our simulations. Some technical complications 
arise in the computation of the second member on 
the right hand side of Eq. (^) . In the evaluation 
of {Sm)sm+Sf^ the use of C*-periodic boundary 
conditions in time direction for the gauge fields 
requires C* boundary conditions in temporal di- 
rection also for fermionic variables (in addition 
to the usual antiperiodic ones), in order to en- 
sure gauge invariance of the fermionic determi- 
nant. This implies relevant changes in the formu- 
lation and implementation of the HMC algorithm 
which are explained in detail in Ref. . 

We have chosen the Polyakov line as the local 
adjoint operator which defines the abelian pro- 
jection. Actually, calling L{x,t) the Polyakov 
line starting at point {x, t), the abelian projection 
is defined by the operator L{x,t)L* {x,t), which 
transforms in the adjoint representation when us- 
ing C* boundary conditions. 

The use of a modified gauge action also im- 
plies changes in the molecular dynamics equa- 
tions. One has to maintain constant the modi- 
fied hamiltonian containing Sm ■ A change in any 
temporal link indeed induces a change in L(x,t) 
and hence in the abelian projection defining the 
monopole field. Therefore the dependence of Sm 
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Figure 1. Results for p on the 16"^ x 4 (open cir- 
cles) and the 32'^ x 4 (filled circles) lattice. 
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Figure 2. Chiral condensate (open circles) and p 
(filled circles) on the 32^^ x 4 lattice. 



on temporal links is non trivial and the equa- 
tions of motion for the temporal momenta become 
more complicated. Details of the changes neces- 
sary in the molecular dynamics equations will be 
published in a forthcoming paper We have 
used a temporal lattice extent Nt = A and spa- 
tial extents Ng — 16, 32. We have chosen to vary 
the temperature, T — 1/ {Nta{l3,mg)), moving in 
the {P,mq) plane while keeping a fixed value of 
m^/mp. To do this and to extract the physical 
scale we have used fits to the nip and masses 
published in [0. In particular we present here 
results obtained at niT^/mp ~ 0.505. 

4. Results and Outlook 

Results obtained for p on the 16'^ x 4 and 32"^ x 4 
lattices are shown in Fig. 1. A clear peak is 
present in both cases and, as can be seen in Fig. 2, 
it is located exactly at the chiral phase transition, 
at a critical temperature Tc ^ 150 MeV. More- 
over, as is visible in Fig. 1, the peak scales roughly 
by a factor of eight when going from Ng — 16 to 
Ns = 32, i.e. cx N^. Previous studies of fuU QCD 
at a similar set of physical parameters ||ll[| have 
found a crossover rather than a phase transition 
at Tc- From our present results we can draw the 
following conclusions: i) the disorder parameter 
shows a behaviour very similar to that observed 
in pure Yang-Mills theories, thus indicating the 
presence of a transition from a dual supercon- 



ducting to an ordinary vacuum state also for full 
QCD. This transition clearly happens in corre- 
spondence of the chiral phase transition; ii) our 
results seem to scale roughly with a critical index 
V ^ 1/3. Of course we still need to collect more 
data near the critical point and at different values 
of Ng in order to draw more precise conclusions 
about the critical behaviour of p at Tc- 
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